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Figure 1: The process 7|,(gi) 'jlil^) Pl(^i) Pl(^2)- The box represents the BFKL 
Green's function, and the t-channel gluons couple to the quark hues in all possible ways. 

1 Introduction 

The next generation of e~^e~ colliders will offer the possibility of clean testing of QCD 
dynamics. By selecting events in which two vector mesons are produced with a large 
rapidity gap in between, through the scattering of two highly virtual photons, one accesses 
a kinematical regime in which a perturbative approach is justified. If additionally one 
selects events with comparable photon virtualities, the perturbative Regge dynamics of 
QCD should dominate and allow the use of resummation techniques of the BFKL [1] type. 
In this paper we study these effects in the case of the reaction (see Fig. [1]) 

lliqi) ll{q2) ^ Pdh) pLih) (1) 

in the forward region, i.e. —t = —{qi — = —tmin, where the virtualities Ql = —qf, 
Ql = — g| of the scattered photons play the role of the hard scales. 

This is an interesting process to study, first because the experimental signature is 
simple: double tagged events in high energy e~^e~ collisions [2], with two longitudinally 
polarized rho mesons in the final state. The non-zero virtualities of the photons allow 
studying also the forward region. Second, the theoretical calculation can be performed 
analytically and gives a simple answer which makes it easy to investigate BFKL effects. 
This is partly because in the forward region we are studying, the amplitude for transversely 
polarized photons vanishes identically [3]. 

No data exist yet for the process ([1]) at sufficiently high energies, but some experimental 
data exist for Q2 small [4] which may be analyzed [5] in terms of generalized distribution 
amplitudes [6]. Related reactions have been calculated, both at the Born level and using 
BFKL resummation. The similar process with real incoming photons, 77 — > V1V2 and 
large t where the hard scale is provided by the momentum transfer instead of the photon 
virtualities, was studied in [7,8], and the process •yp VX at large t, has been studied 
both for heavy [9] and light [10] vector mesons. The advantage of the latter process is that 
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there now exists high quahty data from HERA [11], and the BFKL calculations reproduce 
the measured differential cross sections for longitudinal p meson production and for J /if) 
production very well, while the Born level calculations do not work so well [9-11]. At the 
same time the cross section for transversely polarized p mesons is not well understood [10]. 

In a recent paper [3], the scattering amplitude of the process ([1]) was calculated in 
the Born approximation and the feasibility of a dedicated experiment was proved. The 
production amplitude decreases dramatically with t so that its magnitude at t = tmin 
dictates the rate of the reaction. In this paper, we report on a calculation of BFKL effects 
in the leading order approximation at t = tmin- We find that the leading logarithmic (LL) 
BFKL resummation greatly enhances the cross section, as previously observed in other 
reactions [8,12-14]. In addition, we estimate next-to-leading logarithmic (NLL) effects by 
the use of the BLM scale fixing prescription [15] supplemented by the resummation scheme 
of Refs. [16-18]. We show that there is still a substantial increase in the cross section 
compared to fixed-order calculations, but smaller than the LL calculation. It would be 
interesting to compare our estimates with the ones based on the recently available NLL 
impact factor for the process 7* — > of Ref. [19] and the NLL BFKL kernel [20]. 



2 Leading order BFKL forward amplitude 
2.1 General expression 

In the BFKL framework, the amplitude of the process ([1]) can be expressed through the 
inverse Mellin transform with respect to the squared center-of-mass energy s as 

A{s,t)=^s I ^e-^/.(r^), (2) 

where t — tmin ~ —r^, is considered as Euclidean, as is any two-dimensional vector in 
the following), and Y is the rapidity variable, Y = ln(s/so). The minimum momentum 
transfer is here g iven by tmin ^ Q'lQ^/ ^ ['^]* the particular case where — 0, the 
BFKL Green's function can be easily obtained in momentum space [21]. In this case, the 
impact representation for /a;(0) reads 

1 r dk^ dk'"^ r°° 1 / k'^\''' 

^"(°'4/ffr<(W)/^<'^^(^) . (3) 

where the integration over angles has been performed. The functions are the impact 
factors describing the coupling of the BFKL pomeron to vertex 1 or 2, and a, b are color 
indices. The function (x'(z/) is the BFKL characteristic function which is defined by [1] 

uj{iy) = asx{j^) , (4) 

with (is = Q-sNc/'K and 

X{y) = 2^(1) -ilj(^- + iv\-^{^~ iu\ , ^{x) = V\x)/V{x). (5) 
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The general solution for arbitrary values of t is more involved [22] and will not be discussed 
here. We only want to stress that it has a simple expression at the Born level, 



A = is 



(6) 



which has recently been evaluated in Ref . [3] . The special case = of Eq. ([6]) can be 
readily obtained from Eq. ([3]). This relation will be studied in Section [2^21 The impact 
factor Yl PL is given at t = tmin by [3] 



$f (fc, r = 0) = r'Qi I dz zz (j){z 



ab / 



(7) 



where mf = Qfzz (with z = 1 — z), 4){z) is the distribution amplitude of the meson which 
is here given by its asymptotic form (p(z) = 6zz, ancj^l 



C 



ab 



(2vr) 



e S"'' 
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where fp is the p"— meson decay constant. 
Inserting ([7]) and ([8]) in ([3]) we get 



/^(r2 = 0)=47r«2^ 



em 



c 



fpQiQ2 / dzi zizi (j){zi) / ^2:2 2:2^2 0(2:2) 



(9) 



du 



dk^ 



1 



1 



2iu 



dk 



12 



k'^ V rn? k'"^ + w? 



k 



i—2iv 



Let us denote the expressions in the square brackets as /(z, z/) and I{z, —v) respectively. 
It is straightforward to put the propagators on a common denominator and show that 

3^,„ 3 



[zz 



(10) 



Next, the z integrals can also be done. We have 

"1 _3,. 

dzzz^(z)I(z,u) = -{Ql) s+'^'r (f -iz/) r (-| + / dz6[zz)2- 

r (2 + iu) 



-6v^2-2-2- (g?)"^ r (I - tu) r (-i + ^u) 



since the integral over z is just the definition of the Euler beta function. Finally this 
yields the Mellin transform 



f^{r' = 0) =97rXa 



flQiQ2 / dv 



2\-'i/2+iu /^2\-3/2-i!^ 



UJ — Uj{l') 



r^(|-z.)r^ (| + z.)r(-i-..)r(-| + z.) 

r (2 - iu) r (2 + w) 



'121 



^The first factor l/(27r)^ differs from [3] and comes from tfie convention for impact factors. 
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which immediately leads to the final result for the amplitude ([2]) 



TV^ - 1 1 /"°° /O^X" 

(l-z.)r^ (| + z.)r(-i-z.)r(-i + ^.) 

r (2 - iu) r (2 + iz/) 

We define R = Q1/Q2 and write this as 



J(Y,R)= du e'^^"^^ R^''' — 5^ ^ ^ \ V ^ — ^ 



(13) 



^(3,t™„,gi,g2) = ^5^^r7^ — ^— , (m) 

with 



(15) 



We now want to evaluate the integral J. We have three possibilities at hand: (i) 
numerically, (ii) saddle point approximation and (iii) sum over residues of poles of the 
integrand. The first two methods are used to obtain BFKL results in the rest of this 
paper. 



2.2 Born limit 

A very valuable and non-trivial check on the BFKL result as well as the Born result is 
the fact that in the limit Y (no evolution) or alternatively — > (the gluons in 
the ladder do not couple to each other) the BFKL amplitude must reduce to the Born 
level result. In Appendix |X] we show, using method (iii) for the evaluation of the integral, 
that our BFKL result f|T3l) does indeed reduce to the correct Born level results derived in 
Ref. [3], first for the special case R = 1 and the DGLAP-like limit R^ 1, and finally for 
the general case i? 7^ 1. Note that the calculations presented in the Appendix are done in 
a completely different way than in the Born level calculation in Ref. [3], so the agreement 
is very convincing. 

The main results from performing the Born limit are presented in Eqs. ( ISTl) . ( |55l) and 

(EHD- 



2.3 Saddle point approximation of the BFKL amplitude 

For R = 1 the integrand has a saddle point at z/ = 0, but the product of F-functions that 
multiplies the exponential is not very broad. This means the approximation will yield a 
too large answer. But let us try anyway. The expression multiplying the exponential in 
the integrand at u = takes the value vr^/4, so for the case Qi = Q2 we get 



j(F,i) — } — , y>i, (16) 

4 ^14a.C(3)r' 



5 



so the BFKL amplitude is 



A(s,t = = 0. = 0) ~ ^sn"V^^^^^^==. (17) 

However, we can do better than this. We can keep the general case Qi ^ Q2 with 
R = Q1/Q2. The integral is then 



/oo 
-00 



(1^ 



where g {'J = ^ + ii^) was defined in Eq. Expanding the exponent to second order 

we see that the saddle point is shifted: 

uj"(0)Y 

u{u)Y + 2iiy In ~ t<;(0)F + 2i\nRiy + — -^^^ 

uj"{0)Y f 2t\nRY 2\n^R 
h/ + — -— + — -—+cu(0)Y, (19) 

and we can shift the integration variable accordingly to get a Gaussian integral. The 
result is now 

J(y,i?)^ ; ^ exp - ^/^^"^^^J , r»l, (20) 

^ ^ 4 ^14«,C(3)r V 14a,C(3)r;' ^ ' 

and 

, , n • 5 ^ 9(iV,^-l)«.Wp^ e^^'-^^-^ / In^R \ 

(21) 

These approximations can be compared to numerical evaluations of the integral for 
each value of Y . This is illustrated in Fig. [2l We see that the saddle point approximation 
gives the correct asymptotic behavior for large Y and is less accurate for smaller Y. For 
R = 1 the approximate answer is about 40% too big for F ~ 2 and about 10% too big for 
Y ~ 10. For i? = 10 the approximation is slightly better. Note that the Linear Collider is 
likely to test regions with 1^ ^ 5, where the saddle point approximation works relatively 
well. 



2.4 Leading order results 

To obtain definite predictions we need to define our parameters. The p meson decay 
constant and the electromagnetic coupling take the values fp = 216 MeV [23] and = 
1/137, respectively. There are furthermore three parameters in the calculation: in the 
prefactor, which gives the strength of the coupling of the pomeron to the impact factor; 
Os in the BFKL exponent, which gives the strength of the coupling of the gluons inside 
the pomeron; and the energy scale of the rapidity Y. 
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Figure 2: Exact numerical result (solid line) and saddle point approximation (dashed) to 
the integral J(Y, R) for = 1 (left) and R = 10 (right). 

For all cases with a running strong coupling we use a three- loop running ^^(/i^) [24] 
with /x^ = CaQiQ2- Unless otherwise stated, we use this running coupling with Cq, = 1 in 
the prefactor of the amplitude. 

At LL accuracy, is a fixed parameter, i.e., it does not run with the gluon momenta 
in the BFKL ladder. We choose to, however, let it depend on the given Qi and Q2, 
which are external to the pomeron but provide a reasonable choice; we thus choose a.^ = 
—C(s{.QiQ2)- The pomeron intercept is determined by a^, and it is known to be too large 
when comparing to HERA data. Our chosen values give quite large pomeron intercepts, 
but we do not want to artificially suppress the growth by choosing very small values of 
Instead we will see in Section |3] that the growth becomes slower when NLL corrections to 
the BFKL evolution are included. 

The rapidity is defined as 

^ ^ i^^ok) • <^^' 

where cy is a constant that is not constrained at LL accuracy. As discussed in [14] this 
constant is related to the average attained values of 2:1^2 in the process. The authors 
of [14] chose a very small value cy = 0.01. We estimate the corresponding effect more 
conservatively and choose cy = 0.3 for the cross section predictions shown below (see 
also [13]). 

We will now investigate the sensitivity to these choices. Note that all the results shown 
in this Section have been obtained by numerical evaluation of the integral over u and not 
by the saddle point approximation. 

We begin by comparing the energy dependence of the BFKL cross section and the 
Born cross section. Fig. [3] shows the ratio of the differential cross sections da / dt\t=t^^„ 
calculated from BFKL and at the Born level, as a function of the rapidity Y for three 
choices of Qi and Q2, and for a fixed value ctg = 0.2. There is clearly a strong Y- 
dependence, as seen from (pTi) . Note that all of the F-dependence of this ratio comes 
from the BFKL amplitude, since the Born level result is independent of the energy. Fig. H] 
shows the same plot but using our "standard" choice of cxg = ^0!s{QiQ2)- The difference 
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|A(BFKL)|^ / |A(Born)|^ 



100000. : 




2 4 6 8 10 



Figure 3: Ratio of the leading order BFKL cross section ^ to the Born level cross 
section as a function of the rapidity F, using ots = 0.2. The solid curve is for R = Q1/Q2 = 
1, the dashed curve is for i? = 5, and the dotted curve is for R = 50. 
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Figure 5: Ratio of the leading order BFKL cross section ^1 to the Born level cross 

etc I t'tjiifi 

section as a function of the center-of-mass energy a/s, for Qi = Q2 = Q = 2 GeV (left) 
and Qi = 10 GeV and Q2 = 2 GeV (right). The different curves in each plot correspond to 
three different definitions of the rapidity variable, Y = ln(cys/(QiQ2)- The solid curves 
are for cy = 1, dashed curves are for cy = 1/2, and dotted curves are for cy = 2. The 
scale of as is given by = 1. 

between Figs. [3] and H] serves to illustrate the sensitivity to the choice of a^, and shows 
that using a Q-dependent coupling for the pomeron decreases the growth with energy for 
increasing virtualities. 

In Fig. [5] we show the same ratios as a function of the center-of-mass energy ^/s for 

= g2 = 2 GeV in the left plot and for Qi = 10 GeV, Q2 = 2 GeV in the right plot. 
The three different curves represent different choices of the energy scale in the definition 
of Y, corresponding to three different values of the parameter cy. This freedom to change 
the scale introduces an additional uncertainty in the results. In Fig. [6] we show the same 
kind of plot, but varying instead the parameter Cq in the argument of to highlight the 
uncertainty coming from the choice of scale in as- Note that the parameters cy and Ca 
both affect the argument of the BFKL exponential, and thus the energy evolution. 

Finally, in Fig. [7|we show the ratio as a function of Q = Qi = Q2- Thus we see that the 
BFKL prediction differs from the Born level prediction in all kinematical variables, which 
allows testing BFKL dynamics experimentally, and possibly fitting the free parameters. 

We have now seen that there are several uncertainties in the calculated amplitude. 
This has to be kept in mind when viewing the cross section predictions that will follow. 

In Fig. [S] we show the differential cross section da / dt\t=t^-^ as a function of the photon 
virtuality Q in the symmetric case R = 1, i.e. Q = Qi = Q2, for three different energies 
^/s, and in Fig. [9] we show the same cross section as a function of ^/s for three different 
virtualities Q. 

These predictions are made with the parameter choices discussed above and contain 
the corresponding inherent uncertainties. To get some idea of the possible variation in 
the magnitude of the cross section because of the parameters, we plot in Fig. [10] the 
cross section for Q = 2 GeV for the standard parameter choices, and for two extreme 
versions, one where we choose new parameters = 1/2 Ca, c'y = 2cy and one with 
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Figure 6: Same as Fig. [5], but the different curves in each plot correspond instead to 
three different scale choices in a^. The solid curves are for = 1, dashed curves are for 
Ca = 1/2, and dotted curves are for Cq, = 2. The scale of Y is given by the standard 

Cy = 0.3. 




Figure 7: Ratio of the leading order BFKL cross section ^1 to the Born level cross 

etc I tjYiin 

section as a function of the virtuality Qi for fixed Q2 = 2 GeV, for three different energies 
^/s and standard scale choices Ca, cy as defined in the text. 



10 



dcr/dtl,.,^, [pb/GeV^] 




Figure 8: Differential cross section da / dt\t=tr,un ^ function of tlie plioton virtuality 
Q = Qi = Q2 for center-of-mass energies ^/s = 100 GeV, 300 GeV and 500 GeV, with 
standard scale choices as defined in the text. 




Figure 9: Differential cross section da /dt\ function of the photon-photon centcr- 

of-mass energy y/s for photon virtualities Q = Qi = Q2 = 2 GeV, 3 GeV and 4 GeV. 
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Figure 10: Differential cross section da/dt\t=t^i„ for photon virtuality Q = Qi — 
Q2 — 2 GeV for standard parameters Ca, cy (black curve) and for new parameters 
= 1/2 Cq, c'y — 2cy and c'^ — 2ca, (^y — 1/2 cy (upper resp. lower gray curves). 



c'^ = 2ca, c'y = 1/2 Cy. Tliese curves are plotted in gray and should give some indication 
of the theoretical uncertainty. 

3 Estimation of next to leading order effects 

The BFKL kernel is known to NLL accuracy [20], and the NLL impact factor for our 
process has recently been computed [19]. The calculations to obtain the full cross section 
are difficult, however, and have not yet been performed. It will therefore be interesting to 
estimate the NLL effects on our calculation. When the cross section is finally computed 
it may very well be the first complete NLL cross section to be obtained, and it will be 
very interesting to compare with the estimation presented here. 

To estimate higher order effects we thus implement two improvements to the LL BFKL 
amplitude. First, we use ELM scale fixing [15] for the running of the coupling in the 
prefactor. Second, we use a renormalization group resummed BFKL kernel, as will be 
explained below. 

In Ref. [14] it is shown that for the BFKL calculation of the total 7*7* cross section, 
the BLM procedure for choosing the scale leads to /x^ = CaQiQ^ with Cq = e~^/^, and 
thus to a larger coupling which will increase the cross section. In the process 7*p Vp 
at ncxt-to- leading order [25] the correct scale choice was instead found to be Cq, = e~^/^. 
The use of the BLM procedure in exclusive processes has been further discussed in [26]. 

An approximate BLM scale for our process is found by using the NLL impact factors 
computed in [19] and neglecting higher order effects in the BFKL kernel. The BLM 
procedure then, choosing the scale such that the terms proportional to /3o vanish, leads 
to the simple choice Cq, = 1. 

The NLL kernel, taken as it is, is larger than the leading order kernel, and for any 
reasonable value of ag it is negative; the pomeron intercept becomes less than one for 
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Figure 11: The characteristic BFKL functions ujnll{i) for the NLL resummed model 
with n/ = 3 (sohd hne) and rij = (dotted hne) and 0^(7) for LL BFKL (dashed hne), 
for a value of cxg = 0.2. The dot shows the location of the saddle point of oJNLhiP'f = 3), 
7s = 0.51, while ^(7) and uj^Lbinf = 0) have saddle points in 7 = 1/2. 



as ^ 0.15. Furthermore it has two complex conjugate saddle points which can lead 
to oscillating cross sections. The perturbative expansion of the BFKL kernel is therefore 
highly unstable and far from converging. See [16, 17,27] for a discussion of these problems. 

Different groups have proposed curing the problems of the NLL kernel by certain 
resummations and methods to stabilize the expansion. In particular some resummation 
schemes have been proposed [16-18,28] which lead to modified BFKL kernel^ xil^^) 
depending on both 7 and u. This form comes about because of a resummation which 
removes unphysical double logarithms in the DGLAP and "anti-DGLAP" limits where 7 
is close to or 1 [16], and which resums logarithms from the running of the coupling [17]. 
This means that when performing the inverse Mellin transform in Eq. ([2]), the position of 
the pole in the u plane is determined by the equation 

^ = asXil,^)- (23) 

This equation implicitly defines the function uji\ill{i) that we need to perform the integral 
J{Y,R) over u = i{l/2 — 7). We therefore estimate the NLL BFKL result by using a 
resummed BFKL kernel in the BFKL exponential, but we keep the LL form of the solution 
and the impact factors. 

The renormalization group (RG) resummed BFKL kernel of Salam [16] and Ciafaloni 
et al. [17] is actually a resummation of the full NLL kernel, which removes the problems of 
this kernel. It is however possible to perform such a resummation of the LL kernel, which 
leads to a result that includes a large part of the corrections coming from the NLL kernel. 

^The discussion is most conveniently performed using the variable j — ^ + iv rather than v. 
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One such model, based on the general idea of the approach of [16,17] (more specifically 
Scheme 4 of [16]), was recently proposed by Khoze et al. [18]. This approach uses a fixed 
strong coupling in the BFKL kernel which is essential for the approach here. A running 
coupling in the BFKL kernel radically changes the properties of the solutions, and it is 
no longer possible to evaluate the uj integral by a simple residue. Therefore we do not 
pursue it here, although it may be attempted along the lines of [29]. 

We make one modification to the approach of [18]: they use an asymmetric scale choice 
in the definition of the rapidity (see [16, 17] for a discussion) which was appropriate for 
their problem under study, but our process is more suited for a symmetric scale choice 
and we therefore perform the necessary modification. The characteristic function is then 
expressed as 

X(7,^) = Xo(7) + asXi(7,^) (24) 

where Xo(7) = xil) is the usual BFKL function ([5]). The correction piece xi(7,ci;) is 
given by 

where Xo* is the higher twist part of xo, 

Xo*(7) = Xo(7) - ^ - = 2^(1) - ^(1 + 7) - ^(2 - 7)- (26) 

Ai{uj) is obtained from the Mellin transform of the DGLAP splitting function Pgg by 

^-P,,(^) = - + Ai(^) (27) 



with 

AJu) = - 

12 18 V36 6 

We will in the following throw away any terms proportional to uj and only keep singular 
and constant terms. For = we then have Ai{uj) ~ —11/12. It is possible to account 
for quark loops for ny > by replacing [18] 

r a 1 1 




(29) 

The characteristic kernel UNLiil) obtained by solving fl23l) with x{l,^^) given by fl24l) 
is shown both for nj = 3 and riy = in Fig. [11] together with the LL BFKL kernel. An 
important feature of the kernel fl23l) is that for rij = it has no pole at 7 = 0. The 
pole reappears when including the quark loops as shown above. It is also clear that for 
the resummed kernel both the pomeron intercept and the second derivative at the saddle 
point are reduced. 
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To compute the cross section using this model we need to make use of Eq. (1231) in 
performing the integral over v. This is possible to do purely numerically, by solving the 
equation (1231) explicitly for each given value of 7 when performing the integral numerically. 
However, we may obtain some more insight into the properties of the NLL corrections by 
instead performing the integral by the saddle point method as in Section I2.3[ We expect 
the accuracy of the saddle point method to be similar to the LL calculation. For the 
saddle point calculation we need only compute the position of the saddle point (7s,a;s), 
where Us = uiNLiils), and the particular value u'^ = uj'^l^^^'Js)- These can be obtained 
explicitly in terms of the partial derivatives of the function xil,^) as follows. 

The chain rule gives 



d'y \ 97 duj d'y 



(30) 



so that the saddle point condition gives 

duJNLLil) 



d'y 







1 = 13 



d'y 



0. 



(31) 



7=7s 



We then obtain the second derivative of ujnll{,i) (using the saddle point condition) 



„ d'^UNLhin) 



dY' 



l=ls 



7=7. 1 - 



(32) 



Finally, to find the saddle point values 75,0;^ we must simultaneously solve the equations 



7=7. 
= asX{ls,^^s)- 



0, 



(33) 



The symmetric scale choice with -Uf = gives a kernel symmetric under 7 ^ 1 — 7, so 
that the saddle point is always located at 7s = 1/2. This does not hold for an asymmetric 
kernel obtained either by the asymmetric scale choice, or, as here, by the inclusion of 
n/ > effects. The saddle point values 7s,(x's must then be determined separately for 
each value of a^. 

The residue obtained by the integral over u using the a;-dependent kernel is different 
from the LL case; we have 



Res z — 7 r 



1 - a;,. 



where we defined the third constant 



a;, = a. 



duo 



(34) 



(35) 
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Figure 12: Saddle point approximation of the integral J(Y, R 
the NLL corrected kernel. 



1) for LL BFKL and for 



The saddle point evaluation of the integral JNihiX^R) now gives 



A^^^(s, t — tmini Qli Q2) 



IS vr 



^27r 



9(iVf - 1) e-=^ 



2 In^ R 



m QlQl 



1 - cu.. 



(36) 



We show the difference between the energy evolution of the LL kernel and the re- 
summed NLL kernel in Fig. [121 This plot shows the integral J{Y, R = 1) for LL and 
NLL for a fixed value as = 0.2. The growth with rapidity is strongly reduced by the 
NLL effects, and the diffusion pattern is also changed. This is quantified by the pomeron 
intercept ttp = Ug, which is reduced from ap = 0.55 to ap = 0.20, and by the second 
derivative u'^ of the kernel, which decreases from u'^ = 28asC(3) ~ 6.73 at the LL level to 
uj" ~ 1.02 using the NLL approximationjf] The overall normalization is also affected by 
the factor 1/(1 — Ug), where in this case uJs — —1.51. 

The complete cross section prediction from the NLL modified amplitudes with BLM 
scale choice is shown in Fig. [13] as a function of the energy. The dashed line shows the 
LL BFKL result for comparison. 



4 Conclusions 

Since we expect the linear collider to be able to cover a quite large region in rapidity, exper- 
iments will allow us to test the dynamics of Pomeron exchange. The Born approximation 
estimate [3] of this reaction showed that the process 7*7* pp should be measurable by 



■^Note that these values of uig and cj" 
them correctly when using formula (|36p 



depend on the value of as , so care has to be taken to determine 
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dcr/dt|t=t.^„ [pb/GeV^] 
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Figure 13: Cross section for LL BFKL (dashed lines) and for the NLL corrected kernel 
(solid lines), using cy = 0.3 and the BLM scale choice Cq, = 1, for the three cases Q = 
Qi = Q2 = 2 GeV, 3 GeV and 4 GeV (from top to bottom in the plot). 

dedicated experiments at the next linear collider, for virtualities of the photons up to a 
few GeV^. The increase of the amplitude that we have obtained when considering BFKL 
resummation effects is good news for the experimental feasibility of this study. At leading 
order the growth with energy is very large, but the next to leading order effects seem to 
moderate the growth. 

This exclusive diffractive reaction may as anticipated become the best tool to inves- 
tigate the perturbative picture of the hard Pomeron. For example, considering that the 
full next-to- leading order BFKL calculation of the cross section is likely to appear in 
the near future [19], measurements of this process would be able to test the impact of 
the NLL corrections and the various approaches proposed in the literature. Finally, it 
would also be interesting to pursue this line of research with Odderon exchange processes 
such as 7*7* — * vr^vr'' or, through interference effects in charge asymmetric observables in 
7*7* 7r+7r"7r+7r" [30] 

Note added: Shortly after submitting this paper, the paper [31] appeared, containing 
the full NLL calculation alluded to above, and confirming our LL results. We also note 
that the BLM scale we obtained under our assumptions in Section [3] is not changed if the 
NLL effects on the BFKL kernel are included. 
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A Born limit 



In order to simplify the presentation that follows, wc will make this check in two steps. 
First, we evaluate the integral J(0, 1) exactly, that is, in the case that i? = 1, by picking 
up all the residues of the integrand in the correct half-plane. Second, the method is 
generalized to the case where R is arbitrary. In the asymptotic limit 1, the leading 
contribution can then be easily extracted, as was also done in [3], and further, an explicit 
twist expansion can be written. 

For fixed R, the integral is now just a constant number. Changing to the integration 
variable 7 = | + ii/, we get 



1 /-i/a+ioo T 

^ Jl/2-ioo * 



where 



with 



g{^) = R^^ g,{^) 
7r^7(l — 7) 



Hi) 



r(5/2 - 7)r(3/2 + 7) sin^(7r7) sin3(7r7) ' 
In the following, we will use the auxiliary functions / and /i defined by 

,(^)^_^^i,2._Ai7L_ 
sm (7r7) sm (7r7) 



(37) 

(38) 

(39) 

(40) 



The function g{'^) obviously has triple poles at all integers as given by the sin^(7r7) factor, 
except at 7 = 0, 1 where it has double poles. 

Let us start with the case R — 1. The integration contour (in 7) is parallel to the 
imaginary axis and crosses the real axis between and 1. The contour can be closed 
arbitrarily to the right or to the left. We choose to close the contour in the left half-plane 
and thus to pick up all the poles at 7 = — n, n = 0, 1, 2, . . . 

To extract the triple pole residues at7 = — n,n = l,2,... we write ^ — x — n where 
n e Z, n 7^ 0, 1 and x e R, x <S 1 so that 



f{x - n) 



sm TT x 



n] 



[Tlx) 



1 + 



f{x - n) 



(41) 



The residue at 7 = —n is the residue at a; = 0, which is by definition for the triple poles 

1 (f ( J-11 



Res 0(7) = lim — -— r \ x'" , 



1 + 



f{x - n) 



(42) 
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Taylor expanding f{x — n) around x — Q now easily yields the needed expression 

Res ^(7) = til! (/"(_^) + vrV(-n)) . (43) 

7=-n ZT^^ 

This is the first step in the evaluation. We also need a convenient form for f"{—n) 
that can be summed over n. In order to proceed we rewrite /i without any F-functions 
as 

^^^^^ (l/2-7)(l/2 + 7)(3/2-7) ^ ^ 

Prom this expression, we get after some algebra the derivatives of /i up to the second 
order, at the pole — n, in the form 

n_ . _ , ,,„7r' [ 3 2 3 ^^ 

- ^ ^) 2 lv(2n-l)2 + (2n+l)2 + (2n + 3)2j' 

= -V.(-n) + (-ir2.^(^ + ^ + ^), 
which immediately leads, for n > 0, to 

3 



Res 5(1(7) = - 

7i=-n TT 



+ TT^T + 



.(2n-l)3 ' (2n + l)3 ' (2n + 3)3J ' ^^^^ 

The sum over all residues in the left half plane is therefore related to the definition of the 
Riemann C, function, so 

^ „ , , 63C(3) -46 

Res 51(7) = . (47) 

^-^ 7=-n yTT 
n=l 

To evaluate the integral we should also include the double pole of 5(1(7) in 7 = 0. Defining 

^1(7) = 7' 51 (7) (48) 

one gets 

/ii(0) = |- and ^'i(0) = -^ (49) 
and thus the residue is — 8/(97r). So finally, 

^ Res ^1(7) = ^-^^^ . (50) 

' 7=— n Tf 

n=0 

For our first goal which was to consider = 1, for which 5(7) = 51(7), this sum of residues 
of gi is enough and thus J(0, 1) = 14^(3) — 12. The answer for the Born amplitude at 
Qi = Q2 = Q therefore is 

/V2 _ 1 077-2 

A^°^°(S, tmin, Q) = i'^j^ Ocla^mfl -Q^i^^Ci'^) " 12) (51) 
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which is exactly the same as in [3]. 

Let us now study the general case. Note first that the limits R ^ 1 and i? <^ 1 
are of special physical interest, since they correspond to the kinematics typical for deep 
inelastic scattering on a photon target described through the coUinear approximation, i. e. 
the usual parton model [32]. From Eqs. fl371 - [39|) . the 7-contour can be closed on the left 
(resp. right) half-plane for i? > 1 (resp. R < 1). Let us investigate the case R > 1. 
As above for gi, the integrand g has infinitely many poles: a double pole at 7 = 1 and 
triple poles at all other integers. Since we are closing the contour in the left half-plane we 
need only worry about the poles at 7 = — |n|. Using the method that we used above for 
R = 1 gives an expression for the residues at theses poles, and a related series expansion 
in powers of l/R^ for J{0,R). It consists of a term corresponding to the double pole at 
(which is leading at large R), plus terms that are suppressed with at least a power 
l/i?^" compared to the 7 = pole. Note that this expansion with respect to l/i?^, for 
fixed values of Q2, corresponds to a twist expansion. The full resummation of this twist 
expansion can be performed, as we explain now. 

Denoting 

M7)=/?'^/^i(7) (52) 
one can readily extract the pole of (7(7) at through the relation 



which leads to 



and thus to 



/i'(0) = 21ni?/ii(0) + /i'(0), (53) 



/.'(O)^i^(lnfl-l,, 



J^0,R»1)^'4C^~^.]. (54) 



3 V ^ 6i?y 

which, inserted in formula (|T^ for the amplitude, gives 



A^°™(^, t™n, Ql = RQ2,Q2) = ^^^^ «s«e™/p ^ ' ' ^^^^ 

This is again exactly the result as in [3]. Let us go further. Combining Eqs. P3|) and 
(I38H1QD, one obtains 

Res g{^) = Res g,{^) + {-l)-^^{\n' Rf,{-n) +\nRf[{-n)) (56) 

Using the result (H5|) for /i and f[, and relating the obtained series to the Li2 and Lis 
series 

Li2(^)=Ep' Li3(.) = 5:|„ (57) 

fc=0 k=0 
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1 

R 



one finally gets, after using Eq. (H6|) and ([5 
J(0,i?) = -6{ R 

+ 

+ { Li2 



-6 ( i?- j lni?-3 (r+^ 



In^R 



(5^ 



3 
R 



2 + 3R 
1 



R 



Lis 



ln(l + -J-ln^l 
lni?+ ( Lig 



1 

R 



1 

1 

R 



In^R 



Li, 



1 



which is in agreement with the result obtained in [3], although the formula given there is 
expressed in a different way. After some transformations using the Landen relations [33] 
for Li2 and Lis one finds that the two expressions match. The obtained formula (1581) is 
very convenient to obtain the twist expansion. We illustrate this in the following formula, 
by giving the first twist correction to the amplitude. 



A 



Born / 



5, 



Qi 



IS 



7V2 



RQ2, Q2 

967r2 r 



(59) 



1 
R 



Ini?- 



1 2 „ 47, „ 1307\ 

— - In^ i? H n H + 

I 5 ^75 ^ 2250 / 
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